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Recently Groh and Dietrich QXJ] calculated orientational order transitions 
for dipolar fluids contained in spheroidal shapes of aspect ratio k. Accord- 



ing to their study the thermodynamic state of a dipolar fluid depends on 
Os 



the shape of the fluid's container. For example, a homogeneous fluid in a 

short fat container would phase separate when transferred to a tall skinny 

container of identical volume and temperature. Their calculation thus lacks 
O' 
^ ' a thermodynamic limit. 



Existence of a thermodynamic limit for dipolar fluids has never been 



proven. Indeed, the long range anisotropic dipole interaction presents signif- 
icant difficulties. These include mathematical difficulties such as the condi- 
tionally convergent, shape dependent integrals examined by Groh and Diet- 
rich. Physically, the long interaction range causes domain or texture forma- 
tion to avoid demagnetizing fields. Still, Griffiths [fj] has proven the existence 
of a shape independent thermodynamic limit for the free energy of dipoles 
on lattices. So we presume that dipolar fluids also possess a thermodynamic 



limit. 

The shape dependence in Groh and Dietrich's calculation results from 
their assumption of spatially uniform magnetization. Real crystalline fer- 
romagnets form domains to avoid stray "demagnetizing" fields which arise 
when the magnetization is not parallel to the sample boundary. In a liquid 
ferromagnet, however, the domain wall thickness is limited only by sample 
sizef^J. Thus the magnetization is expected to rotate across the sample in a 
manner as yet undetermined. There is one sample shape which avoids the 
demagnetizing field - the limit k —> oo corresponding to an infinitely pro- 
late spheroid. The bulk free energy calculated for a uniform magnetization 
in this needle-like shape will match the true free energy of any shape with 
appropriate spatially varying magnetization. 

While the k — > oo limit is the uniquely appropriate value for calculating 
bulk free energy density assuming a spatially uniform magnetization, we do 
not expect a real ferromagnetic liquid drop to take this needle-like shape in 
the absence of external field. The actual shape and magnetization texture of 
a ferromagnetic liquid droplet remains an important unsolved problem. 

Taking k — > oo is not practical for computer simulations. Recent com- 
puter simulations [4-7] solve the long range interaction problem by combining 
Ewald summation with a reaction field ||. The limit of infinite surrounding 
dielectric constant e — > oo cancels the demagnetizing field and yields the 
true thermodynamic limit with a uniform state, while being computation- 
ally more tractable than taking k — > oo. In simulations with e = 1, Weis 



and LevesqueQ] find domain formation consistent with our assertion that 
demagnetizing fields create spatially nonuniform states. 

Now consider Groh and Dietrich's density functional theory phase dia- 
gram in the k — >• oo limit. This contains a continuous magnetic transition at 
high temperatures and a tricritical point below which phase separation oc- 
curs between dilute isotropic gas and dense magnetized liquid. Simulations 
of hard@] and soft[j|, |(| dipolar spheres place the magnetized liquid state at 
higher densities and lower temperatures than suggested by the present cal- 
culation. A simulation of the Stockmayer fluid|7] reveals isotropic phase sep- 
aration with a conventional critical point. We believe the density functional 
theory may be brought into closer agreement with computer simulation by 
incorporating effects of random particle positions which shift the magnetic 
transition to larger densities^, and particle chaining||10| which competes 



with ordinary phase separation at strong dipolar coupling. Both effects re- 
quire three-body and higher terms in the density functional theory. 

The phase diagrams found by Groh and Dietrich illustrate part of a 
generic sequence of phase diagrams for dipolar fluids|TT|. Many parame- 
ters, including the strength of the 1/r 6 attraction in the Stockmayer fluid, 
carry the phase diagrams through such a sequence. But sample shape is not 
among such parameters. 

This work was supported by NSF grant DMR-9221596 and by the A. P. 
Sloan foundation. We are indebted to R.B. Griffiths for numerous discus- 
sions. 
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